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Abstract 

A general purpose math model of a servovalve controlled hydraulic actuator system is 
derived. The system consists of a linear actuator with unequal piston areas, a single stage 
servovalve, a gas charged hydraulic accumulator, and the interconnecting piping. The state 
equations are integrated using the Advanced Continuous Simulation Language (ACSL) for 
determining the system’s dynamic response characteristics. Using this generalized hydraulic 
actuator system model, response characteristics were determined for various servovalve 
commands. 

Nomenclature 
Variables 
A area 

B fluid bulk modulus 

C viscous damping 

DP differential pressure 

/ servovalve bandwidth 

Fst static force 

Fstk sticking force 

Fcol coulomb friction force 

K pressure loss coefficient 

m fluid mass 



M 

total mass of piston and load 

P 

pressure 

Q 

volumetric flow 

s 

laplace operator 

t 

time 

V 

volume 

Wn 

frequency 

X 

actuator displacement 

y 

servovalve spool displacement 

7 

ratio specific heats 

T 

servovalve time constant 

p 

fluid density 

( 

damping 

Subscripts 


acc 

accumulator 

act 

actuator 

c 

actuator cap side 

9 

gaseous 

r 

actuator rod side 

cruel 

command 

in 

initial condition 

l 

leakage 

k 

liquid 

rt 

return side 

sp 

supply side 

tot 

total 
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vrt servovalve return side 

vsp servovalve supply side 

v servovalve 

va servovalve port a 

vb servovalve port b 

vc servovalve cap side 

vr servovalve rod side 

1.0 Introduction 

Many dynamic systems use linear or rotary hydraulic actuators to generate large forces 
for obtaining high speed response. To satisfy design requirements for such systems, it is 
necessary to accurately predict load position, velocity, acceleration and system stability. 

A generalized hydraulic actuator system was modeled utilizing the Advanced Continuous 
Simulation Language ( ACSL) computer code to solve the system of nonlinear time-dependent 
ordinary differential equations. Model parameters are easily changed for purposes of 
component and system evaluation. Control strategies for particular hardware components 
may be evaluated by assessing the effect on system performance. For a given operating 
condition and controller gain setting, ACSL evaluates the poles of the system transfer 
function, thus providing a method to assess the stability of the system. 

Three simulations were performed using typical system parameters; a step command 
to the servovalve, a steep ramp command, and a shallow rarrip command. Results of the 
simulations are discussed with selected variables plotted. 

2.0 Governing Equations 

Shown in Figure 1 is a schematic of a typical servovalve controlled hydraulic actuator 
system. The system model consists of an actuator, single stage servovalve, gas charged 
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accumulator, and interconnecting piping. Presented here is the development of the governing 
equations for each component. The model of the hydraulic system presented here is an 
extension of that presented in reference (1) and in addition includes an unequal area piston 
actuator, Coulomb friction, sticking and pressure losses in the piping system. 

2.1 Actuator 

The actuator is modeled as a second order spring-mass system with viscous damping, 
Coulomb friction and sticking. The cap and rod side piston areas are unequal. Piston seal 
leakage is assumed laminar and given by: 


Ql = I<act(Pc-Pr) ( 1 ) 

Referring to Figure 1, continuity of mass flow into the cap end of the actuator yeilds: 

dm c „ _ 

+ QiPl = Qvcp c (2) 

The mass contained in the cap end is given by: 


Differentiating, 


m c = p c V c 


dmc dVc dp c 

nr = Pc nr + Vc nr 


Fluid bulk modulas is defined as: 


Volume rate of change is given as: 


B = 


dP 

dp/p 


dVc dx 

— : — = Ac— 
dt dt 


Substitute (4), (5), and (6) into (2) yeilding: 


Qvc 


VcdPc . , dx 


B dt 1 dt 


i^-n + Qi 


( 3 ) 

( 4 ) 

( 5 ) 

( 6 ) 

( 7 ) 
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Similarly for the rod end: 


( 8 ) 


VrdPr dx 

Qvr = tit ~ Ar Tt ~ Ql 


A force balance on the actuator yields: 


r/" t dor 

M —77 = PcAc - PrAr - C— - (Fst + Fstk + Fcol) 
dt 2 dt 


( 9 ) 


At the point of zero actuator velocity, the fluid pressure must overcome the breakout or 
sticking force, Fstk, before motion is established. At that time a constant force due to 
Coulomb friction, Fcol, opposes the direction of motion. These non-linear terms have the 
effect of increasing system damping. 


2.2 Servovalve 

The servovalve spool actuator is 


modeled as a first order system with response given as: 
T^Ay^Ycmd ( 10 ) 


with 

r = 1/(2 tt/) 

Servovalve bandwidth, f, is typically specified by the manufacturer. The oil flow of the 
servovalve is a function of the square root of differential pressure across the valve and is 
given as; 


Qvc = Kv(P vsp - P vb ) l/ ' 2 

(11) 

<C 

'S 

II 

<s 

c 

1 

(12) 


The variable, Kv, is computed from manufacturer ratings of flow for a given differential 
pressure. For instance, if the servovalve is rated for 40 gpm at a pressure loss of 1000 psi, 
Kv would be calculated as: 

(13) 


Kv = (3.85)(40)/(1000) 1/2 = 4.87m 4 /(sec(/6/) 1/2 ) 
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A servovalve spool may be manufactured with overlap to insure positive shutoff. An 
overlapped spool results in a deadband about the neutral position and acts to increase 
system damping. Deadband is modeled with an ACSL command by setting the dependent 
variable to zero when the independent variable lies between a specified lower and upper 
bound, (ref. 2) 

2.3 Accumulator 

The accumulator is charged with gas which is assumed to behave isentropically. Accu- 
mulator pressure is shown to be: 


Paccin(V gin? = Pacc(Vg) 7 

(14) 

Vg = Vtot - V lq 

(15) 

dV l q n 


a ~ Qsr 

(16) 

^Iq = J Qspdt 

(1?) 

Pace = Paccin{V9in/(V to t ~ [ Qspdt)? 

(18) 


2.4 Piping 

Pressure losses in the inncrconnection piping is modeled by use of pressure loss coefficients. 
There are four significant sections of piping; the accumulator to servovalve, servovalve to 
return, and servovalve to rod and cap ends of the actuator cylinder. These equations are 
defined as: 


DP S p = Ksp(Qs P ) 2 

(19) 

DP C = K c (Qcf 

(20) 

DP r = I<r(Qr) 2 

(21) 

CM 

II 

cf 

Q 

(22) 
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2.5 Transfer Function 

The system of equations for the generalized hydraulic actuator system may be linearized 
by making the following assumptions. 

1. Coulomb friction and sticking are negligible 

2. Pressure losses in the piping system are negligible 

3. Supply pressure remains constant 

4. Zero servovalve overlap 

Incorporating these assumptions, a block diagram of the linearized system is shown in figure 
2. The block diagram may be reduced to yeild the system transfer function: 


A" 


Y cn,d S (t S+1)[S 2 + &S + $($ + t|) 

Which may be put in the form; 

x 

Ycrnd S(tS + 1) [S 2 + 2 (W N S + W 2 ] 

yeilding system natural frequency and damping. 

*2 a2' 


(23) 


(24) 


K 


B_ 

M 


A; A 2 \ 

V r + Vc) 

t = C/(2W N M) 


(25) 

(26) 


For the case A r = A c and V r = V c , such as a double acting linear actuator operating 
about midstroke the transfer function would reduce to: 

WlKvsffylA 


X 

Y cmd ~ S(rS + 1) [S 2 + 2 iW N S + W%\ 

9 2 BA 2 

w » = TW 

( = C/(2W n M ) 


(27) 

(28) 
(29) 
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A rotary motor can also be modeled if the mass of the load is replaced with inertia of 
the rotary motor. 

3.0 Method of Solution 

A computer model of the above equations (Appendix 7.0) was implemented using the 
ACSL code. ACSL was written to facilitate the solution of a system of nth order, nonlinear 
ordinary differential equations which describe a physical system. Special functions are 
available in addition to standard Fortran commands to provide flexibility in the problem 
description. The code sorts the commands into an executable sequence such that a variable 
is calculated before it is used, allowing the system of equations to be inputed in a logical 
sequence. Nine integration algorithms are available (default is fourth order Runge-Kutta) in 
addition to a user supplied routine. The code assigns lower order derivatives a variable name 
then simultaneously solves the system of first order equations and determines the system 
states. 

A post run command invokes a linear analysis which determines the system poles and 
damping. The linearized forcing matrix is determined for perturbations about a particular 
operating point. This capability is particularly useful for designing or developing control 
strategics. 

4 .0 Results 

A typical use for a servovalve controlled hydraulic actuator system would be that of 
model injection into the test chamber of a supersonic wind tunnel. Table I shows system 
parameters for such an application. Dynamic behavior for this model was simulated by 
solving the nonlinear equations described in section 2 using the ACSL code. Results are 
presented for three different servovalve command functions. 
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4.1 Case 1, Step Command to Servovalve 

Figure 3a-e presents results for a step command to the servovalve. The servovalve was 
given a step command of 30%, held for 0.5 second than returned to the neutral position. 
Table I shows the ACSL equation for the step command. Simulated response to this 
command shows large initial acceleration as the servovalve opens rapidly subjecting the 
actuator piston to large differential pressure. While the servovalve spool position is held 
fixed, the actuator piston velocity remains fairly constant, decreasing slightly as accumulator 
pressure decreases. Upon return of the servovalve spool to the neutral position, the actuator 
again undergoes large acceleration/deceleration. Oscillation (second order response) of the 
piston and load occurs due to the compressibility of fluid (fluid bulk modulus). These 
oscillations are damped by viscous and coulomb friction. The maximum flow for this 
simulation occurs to the cap end of the actuator, about -fT5GPM, while the rod end 
experiences about a — 7GPM flow rate (flow in is taken as positive, flow out is negative). 
The difference in these flow rates is due to the unequal piston area of the rod and cap ends. 


4.2 Case 2, Steep Ramp Command to Servovalve 

Figure 4a-e presents results for a steep ramp command to the servovalve. The servovalve 
was given a command of 150%/sec for 0.2 second, then held at 30% for 0.2 second and 
returned to the neutral position at a rate of —150%/sec. Table I shows the ACSL equation 
for tli is ramp command. From these results it is seen that the actuator piston and load is 
subjected to much lower peak acceleration than in Case l. Maximum velocity and flows are 
similar due to the same command hold state. Piston oscillation is much less pronounced 
and of shorter duration than that of Case 1. 
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4.3 Case 3, Shallow Ramp Command to Servovalve 

Figure 5a-e presents results for a shallow ramp command to the servovalve. The servovalve 
was given a command of 60%/sec for 0.5 second, then held at 30% for 0.1 second, and then 
returned to the neutral position at a rate of —60%/sec. Table I shows the ACSL equation 
for this ramp command. Similar to the results of Case 2, peak acceleration of the piston and 
load is again reduced with the less abrupt servovalve spool movement. Maximum velocity 
and flows are again about the same due to the same command hold state as the previous 
simulations. Piston oscillation is almost non-existent for this case. 

These simulations indicate that the acceleration of the piston and load are determined by 
the rate of change of the servovalve command. The magnitude of the servovalve command 
determine piston and load slew rate. 

5.0 Conclusion 

A general purpose math model of a servovalve controlled hydraulic actuator system is 
derived for use with the Advanced Continuous Simulation Language (ACSL). Parameters 
typical of a hydraulic system were input and the model exercised. The math model is of 
a general hydraulic actuator system which can be easily changed. This allows for quick 
assessment of alternate system hardware or control strategy. 
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Table I Example Simulation Parameters 


Cylinder Length 

Actuator 

17 

Inch 

Stroke 

15 

Inch 

Mass of Piston & Load 

2050 

Lbm 

Cylinder I.D. 

2.0 

Inch 

Rod Diameter 

1.375 

Inch 

Kc[ac< 

0.0 

in 5 /(sec Lbf) 

Fstk 

75 

Lbf 

Fcol 

50 

Lbf 

Bandwidth 

Servo Valve 
25 

Hz 

Max Flow @ 1000 psid 

40 

GPM 

Deadband 

+/- 2 

% 

Accumulator and Piping 


Volume 

5 

Gal 

Initial Gas Charge 

1800 

psia 

Supply Pressure 

3000 

psig 

psig 

Return Pressure 

25 

Kqsup 

0.0125 

Lbf sec^/in® 

Kqret 

0.0125 

Lbf sec^/in® 

Kqcap 

0.0125 

Lbf sec^/in® 

Kqrod 

0.0125 

Lbf sec^/in® 


Command to Servo Valve 


Case 1. (Figure 3), Step Command 
Y cm<t = 30(step(.2)-step(.7)) 

Case 2. (Figure 4), Steep Ramp Command 
Y cmd — 150(rainp(.2)-ramp(.4)-ramp(.6)+ramp(.8)) 
Case 3. (Figure 5), Shallow Ramp Command 
Y cm j = 60(ramp(.l)-ramp(.6)-ramp(.7)+ramp(1.2)) 
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Load 


+ x 
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Figure 1. Schematic of Servovalve Controlled Hydraulic Actuator System. 






Note: P sup = Supply pressure: Return 
pressure assumed zero 


Transfer function: 


X 

Ycmd 



s<ts*i)[s’ + £s ♦£($$)] 


For A r = A c = A 
and V r = V c = V 
Transfer function reduces to: 


X 

Ycmd 


2A BK vy gpT /(VM) 
S(xS+1) (S 2 S +2BA 2 /(VM)) 


Figure 2. Block Diagram of Linearized Servovalve Controlled Hydraulic Actuator System for 
Small Perturbations. 


14 








Figure 3a. Servovalve Command (YCMD) and Actuator Position (X) Versus Time for a Step 
Command to the Servovalve. 


c n 


O 

Q 

X 



Figure 3b. Actuator Velocity (XDOT) and Acceleration (XG) Versus Time for a Step 
Command to the Servovalve. 
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PROD, psig 



Figure 3d. Rod Side Flow (QRODG) and Pressure (PROD) Versus Time for a Step 
Command to the Servovalve. 
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Time, sec 

Figure 3e. Supply Pressure (PACC) and Supply Flow (QSUPG) Versus Time for a Step 
Command to the Servovalve. 
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Figure 4c. Cap Side Flow (QCAPG) and Pressure (PCAP) Versus Time for a Steep Ramp 
Command to the Servovalve. 



Figure 4d. Rod Side Flow (QROD) and Pressure (PROD) Versus Time for a Steep Ramp 
Command to the Servovalve. 
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Figure 4e, 
Ramp 




Figure 5a. Servovalve Command (YCMD) and Actuator Position (X) Versus Time for a 
Shallow Ramp Command to the Servovalve. 



Figure 5b. Actuator Velocity (XDOT) and Acceleration (XG) Versus Time for a Shallow 
Ramp Command to the Servovalve. 
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Figure 5c. Cap Side Flow (QCAPG) and Pressure (PCAP) Versus Time for a Shallow Ramp 
Command to the Servovalve. 



Figure 5d. Rod Side Flow (QRODG) and Pressure (PROD) Versus Time for a Shallow Ramp 
Command to the Servovalve. 
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Figure 5e. Supply Pressure (PACC) and Supply Flow (QSUPG) Versus Time for a Shallow 
Ramp Command to the Servovalve. 
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7.0 Code Listing 


PROGRAM CYLINDER 

II __ __ 

*“* „ II 

SERVOVALVE CONTROLLED HYDRAULIC ACTUATOR SYSTEM ANALYSIS 
" ORIGINAL VERSION BY J . F. WATSON, FENGD , FACS 

" REVISED BY 

" D. A. BABCOCK MAY 1989, FENGD, PSS 

II 

— — — — — — — — — — II 

M 

II 

THIS UPDATED VERSION OF THE ORIGINAL CODE INCLUDES THE EFFECTS OF " 
ACTUATOR STICKING AT ZERO VELOCITIY; ACTUATOR COULOMB FRICTION; 

ACTUATOR LEAKAGE FROM THE HIGH TO LOW PRESSURE PORTS; AND 
SERVOVALVE DEADBAND AS A PERCENTAGE OF OVERLAP: THESE NONLINEAR 
TERMS INCREASE THE SIMULATED SYSTEM DAMPING 


|| ************************** VARIABLE TABLE ****************************" 

II 


II 

X 

LOAD DISPLACEMENT 

INCH 

II 

II 

XDOT 

LOAD VELOCITY 

IN/SEC 

H 

II 

XDDOT 

LOAD ACCELERATION 

IN/SEC**2 

II 

II 

II 

IX 

INITIAL POSITION 

IN 

II 

II 

CYSTK 

ACTUATOR STROKE LENGTH 

IN 

ll 

II 

II 

CYDIA 

ACTUATOR PISTON DIAMETER 

IN 

II 

II 

RODDIA 

ACTUATOR ROD DIAMETER 

IN 

II 

II 

WACTR 

ACTUATOR PISTON AND ROD WEIGHT 

LB 

ll 

II 

II 

WLOAD 

LOAD WEIGHT 

LB 

II 

II 

FCYL 

CYLINDER GENERATED FORCE 

LB 

It 

ll 

II 

KFV 

VISCOUS DAMPING 

LB-SEC/ IN 

II 

II 

FVISC 

VISCOUS FRICTION FORCE 

LB 

ll 

II 

FSTAT 

GRAVITY FORCE OF LOAD AND ACTUATOR 

LB 

>1 

II 

FSTK 

ACTUATOR STICKING FORCE 

LB 

It 

tl 

FCOL 

ACTUATOR COULOMB FRICTION FORCE 

LB 

ll 

II 

II 

MTOTL 

TOTAL MASS OF MOVING ITEMS 

LB-SEC** 2/ I N 

II 

II 

PROD 

ROD SIDE CYLDR PRESSURE 

PSIG 

II 

ll 

II 

PCAP 

CAP SIDE CYLDR PRESSURE 

PSIG 

II 

II 

DPPROD 

LINE PRESSURE LOSS ROD SIDE 

F'SID 

ll 

II 

II 

DPPCAP 

LINE PRESSURE LOSS CAP SIDE 

PSID 

It 

II 

QROD 

ROD SIDE CYLDR OIL FLOW 

IN**3/SEC 

ll 

II 

It 

QCAP 

CAP SIDE CYLDR OIL FLOW 

I N** 3/SEC 

II 

II 

VROD 

ROD SIDE CYLDR OIL VOLUME 

IN** 3 

II 

II 

VCAP 

CAP SIDE CYLDR OIL VOLUME 

IN** 3 

II 

II 

AROD 

ROD SIDE AREA 

IN**2 

II 

II 

II 

ACAP 

CAP SIDE AREA 

IN** 2 

ll 
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II 

QL 

ACTUATOR LEAKAGE 

IN** 3/SEC 

II 

QSUP 

OIL SUPPLY 

FLON 

IN**3/SEC 

il 

QRET 

OIL RETURN 

FLOW 

IN** 3/SEC 

II 

PACC 

ACCUM SUPPLY PRESSURE 

PSIG 

II 

POILCG 

INITIAL SYSTEM PRESSURE 

PSIG 

II 

PRET 

RETURN PRESSURE 

PSIG 

II 

VOLACC 

ACCUMULATOR VOLUME 

GAL 

II 

II 

PN2CHG 

INITIAL NITROGEN CHARGE 

PSIA 

II 

Y 

SERVOVALVE 

SPOOL DISP 

% 

II 

YCMD 

COMMAND TO 

SERVO VALVE 

% 

II 

YDOT 

SERVOVALVE 

SPOOL VELOCITY 

%/SEC 

II 

TAU 

SERVOVALVE 

TIME CONSTANT 

SEC 

II 

DBAND 

SERVOVALVE 

DEAD BAND 

% 

II 

SERVFL 

SERVOVALVE 

FLOW AT RATED PRESSURE 

GPM 

II 

II 

SERVPR 

SERVOVALVE 

PRESSURE AT RATED FLOW 

PSID 

II 

PSRVP 

SERVOVALVE 

SUPPLY PRESSURE 

PSIG 

II 

PSRVT 

SERVOVALVE 

RETURN PRESSURE 

PSIG 

II 

PSRVA 

SERVOVALVE 

CONTROL PORT PRESSURE 

PSIG 

II 

PSRVE: 

SERVOVALVE 

CONTROL PORT PRESSURE 

PSIG 

It 

DPPSUP 

LINE PRESSURE LOSS SUPPLY SIDE 

PSID 

II 

II 

DPPRET 

LINE PRESSURE LOSS RETURN SIDE 

PSID 

II 

QVR 

SERVOVALVE 

ROD SIDE FLOW 

IN** 3/SEC 

II 

QVC 

SERVOVALVE 

CAP SIDE FLOW 

IN** 3/SEC 

II 

KQSER 

SERVOVALVE 

FLOW COEF 

IN** 4/ (SEC-LB**. 5 

II 

KQSUP 

SUPPLY PIPING LOSS COEF 

LB-SEC**2/IN**8 

II 

KQRET 

RETURN PIPING LOSS COEF 

LB-SEC**2/I N**8 

II 

KQCAP 

CAP PIPING 

LOSS COEF 

LB-SEC**2/IN**8 

It 

KQROD 

ROD PIPING 

LOSS COEF 

LB-SEC** 2/ I N**8 

It 

II 

KQACT 

ACTUATOR LEAKAGE COEF 

IN**5/( SEC-LB) 

II 

B 

FLUID BULK 

MODULUS 

PSI 

II 

GAMM 

IDEAL GAS SPECIFIC HEAT RATIO 

NA 

II 

ATM 

LOCAL ATMOSP PRESSURE 

PSIA 


it 


'! ic‘kikikicicici<^ci<ikikik^-kic'kic M k^k'k^:'kikikik-k CONST ANT S 

ii 

" SERVOVALVE 


CONSTANT SERVFL = 40.0, SERVPR = 1000.0 
CONSTANT DBAND = 0.02, TAU = 0.006 

II 

" ACCUMULATOR AND PIPING 

ll 


CONSTANT PN2CHG = 1800.0, VOLACC = 5.0, GAMM =1.4 

CONSTANT KQSUP = 0.01250, KQRET = 0.01250 

CONSTANT KQCAP = 0.01250, KQROD = 0.01250 

tl 

" FLUID CONSTANTS 

II 

CONSTANT POILCG = 3000.0, PRET = 25.0, B = 100000.0 

II 

.. ACTUATOR 

25 


CONSTANT CYSTK 
CONSTANT WACTR 
CONSTANT FSTK 
CONATANT ALL 
CONATANT AUL 


17.0, 

CYDIA 

= 2.0 , 

RODDIA = 

1.375 

500.0, 

WLOAD 

= 2000.0, 

KFV 

50.0 

75.0, 

FCOL 

= 50.0, 

KQACT = 

0.01 


0.5 ACTUATOR LOWER INTEG LIMIT ■ 

16.5 ACTUATOR UPPER INTEG LIMIT ■ 


II 


OTHER 


LOGICAL STUCK 
CONSTANT TSTRT 
CONSTANT TSTP 
CONSTANT ATM 
CONSTANT PI 
CONSTANT INPGAL 
CONSTANT CNVRT 
CONSTANT ACCG 
CONSTANT LL 
CONSTANT UL 


* “ 
0.0 
1.0 
14.7 
3.141 

231.0 
3.85 
386.4 

5.0 
5000 . 


TRUE IF VELOCITY IS ZERO OTHERWISE FALSE 
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* M IN** 3 PER GALLON “ 

*'■ IN** 3/SEC TO GPM * 

*" GRAVITATIONAL CNST IN/SEC** 2 
*" LOWER PRES INTEGRATION LIMIT 
%" UPPER PRES INTEGRATION LIMIT 


If 

It 


INITIAL SET INITIAL CONDITIONS " 


ACCUMULATOR 


VLN20 

KACCM 


VOLACC * INPGAL * PN2CHG / ( POILCG + ATM 
(. POILCG + ATM ) * ( VLN20 ** GAMM ) 


II 


ACTUATOR 


RESET ( "NOEVAL" ) 

STUCK = XDOT .EQ. 0.0 
FVISC = 0.0 


ACTUATOR STICKS WHEN VELOCITY IS ZERO " 


FSTAT = WACTR + WLOAD 

MTOTL = FSTAT / ACCG 


ACAP 

AROD 


(PI / 4.0) * (CYDIA ** 2) 

ACAP - (PI / 4.0) * (RODDIA ** 2) 


I PROD 
IPCAP 

IX 


(POILCG * ACAP - FSTAT) / (AROD + ACAP) 
POILCG - IPROD 
0.0 


SERVOVALVE 


II 


II 


KQSER = CNVRT * SERVFL / ( SQRT ( SERVPR )) 

END *"OF INITIAL" 

l***********************************^************^^^^^^^^^* 

DYNAMI C 
DERIVATIVE 
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YCMD = 0 .3*(STEP( .2) - STEP(.7)) H ♦"STEP COMMAND" 

YCMD = 1 .5*(RAMP< .2)-RAMP( .4)-RAMP( .6)+RAMP( .8) ) ♦"RAMP COMMAND" 

|| YCMD = 0 . 6* ( RAMP ( . 1 ) -RAMP ( . 6 ) -RAMP ( . 7 ) +RAMP (1.2))" ♦"RAMP COMMAND" 

II 

"********************** CONTROLLER SECTION ************************** 

II 

II 

|| SERVOVALVE CONTROLLER MAY BE ADDED HERE 

II 

II 

II 

kikikiki<ici(ici<:kirkick:^k-kicik-kik'kicic'kk:^c SERVOUALVE ^ic^icic^icic^icikikicikikicicicicic’k^icic-k-kk'kk: M 

n 

ii 

" DYNAMIC SECTION 

II 

II 

YDOT = (YCMD - YV) / TAU 

YV = LIMINT ( YDOT, 0.0, -1.0, +1.0 ) 

Y = DEAD ( -DBAND , DBAND , YV ) 

II 

II 

" FLOW SECTION 

II 

II 

FLONO = 0.0 

FLOW1 = Y*KQSER*( SQRT ( ABS( PSRVP-PSRVB) ) )*SI6N( 1 . , PSRVP-PSRVB) 

FL0W2 = Y*KQSER*(SQRT(ABS(PSRVB-PSRVT)))*SIGN(1. , PSRVB-PSRVT) 

FLOWS = -Y*KQSER*(SQRT(ABS(PSRVA-PSRVT) ) )*SIGN(1 . , PSRVA-PSRVT ) 

FL0W4 = -Y*KQSER*( SQRT ( ABS ( PSRVP-PSRVA) ) )*SI GN( 1 . , PSRVP-PSRVA) 

QVR = FCNSW( Y , FL0W4 , FLOWO , FL0W3) 

QVC = FCNSW( Y , FL0W2 , FLOWO , FLOW1 ) 

II 

II 

'kkk'kkickkkikiakikikik'k'k ACTUATOR FORCES AND POSITION k'k'kkjrtc'k’kTkjckk-kkkkkic'k' 1 

n 

ii 

X = LIMINT ( XDOT , IX, ALL, AUL ) 

XDOT = INTEG ( XDDOT , 0.0) 

XDDOT = ( FTOTL + FCF ) / MTOTL 

II 

II 

FTOTL = FCYL - FVISC - FSTAT 

FVISC = KFV * XDOT 

FCYL = PCAP * ACAP - PROD * AROD 
FCF = RSW ( STUCK, -FTOTL, FASF ) 

FI A = RSW ( STUCK, ABS ( FTOTL ) - FSTK , XDOT ) 

SCHEDULE STICK .XZ. FI A 

II 

II 

" ***************** ROD SIDE FLOW AND PRESSURE ******************* ** » 

II 

H 

PROD = LIMINT ( PDROD, IPROD, LL, UL ) 

QROD = AROD * XDOT 

PDROD = ( QVR + QROD + QL ) * ( B / VROD ) 

VROD = AROD * ( CYSTK - X ) 

DPPROD = ( QROD ** 2 ) * KQROD * SIGN ( 1.0, QROD ) 

PSRVA = PROD - DPPROD 

II 

II 

QL = KQACT * ( PCAP - PROD ) ♦" ACTUATOR LEAKAGE " 

ll 

it 

kkkkikkk:kkkkkkkkk*k: CAP SIDE FLON AND PRESSURE kkkkkjcjcjrtckkkkkkkkirkk: “ 

a 

ii 

PCAP = LIMINT ( PDCAP , I PCAP , LL, UL ) 

QCAP = ACAP * XDOT 

27 


PDCAP - ( QVC - QCAP - QL ) * ( B / VCAP ) 

VCAP = ACAP * X 

DPPCAP = ( QCAP ** 2 ) * KQCAP * SIGN ( 1.0, QCAP ) 

PSRUB = PCAP + DPPCAP 

II 

****************** ACCUMULATOR / SUPPLY AND RETURN ***************** 

DELVOL = INTEG ( QSUP, 0.0 ) 

V0LN2 = VLN20 + DELVOL 

PACC = KACCM / ( V0LN2 ** GAMM ) - ATM 

DPPSUP = ( QSUP ** 2 ) * KQSUP * SIGN ( 1.0, QSUP ) 

PSRVP = PACC - DPPSUP 

QSUP = FCNSN ( Y, -QROD, 0.0, QCAP ) 

= ** 2 ) * KQRET * SIGN ( 1.0, QRET ) 

PS RUT = PRET + DPPRET 

QRET = FCNSW ( Y, -QCAP, 0.0, QROD ) 

END $ “OF DERIVATIVE" 


DISCRETE STICK 
PROCEDURAL 

STUCK = .NOT. STUCK .AND, 
RSN (STUCK, 0.0, 
0.0 $" SET VEL 0 
"OF PROCEDURAL" 
"OF DISCRETE" 


FASF 

XDOT 

END * 
END * 


ABS( FTOTL 
SI GN( FCOL , 
IF STICKING 


) . LT . FSTK 
-FTOTL)) 


II 


********************** ENGINEERING VARIABLES *********************** 


CONVERT IN**3 TO GALLON PER MINUTE 
QCAPG = QCAP / CNVRT 

QRODG = QROD / CNVRT 

QSUPG = QSUP / CNVRT 

QRETG = QRET / CNVRT 

" CONVERT IN/SEC**2 TO G 
XG = XDDOT / ACCG 

II 

TERMT (T .GE. TSTP) 

II 


END * "OF DYNAMIC" 

END % "OF PROGRAM" 
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